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Abstract. Unstable periodic orbits are known to originate scars on some eigenfunctions of classically chaotic 
systems through recurrences causing that some part of an initial distribution of quantum probability in 
its vicinity returns periodically close to the initial point. In the energy domain, these recurrences are 
seen to accumulate quantum density along the orbit by a constructive interference mechanism when the 
appropriate quantization (on the action of the scarring orbit) is fulfilled. Other quantized phase space 
circuits, such as those defined by homoclinic tori, are also important in the coherent transport of quantum 
density in chaotic systems. The relationship of this secondary quantum transport mechanism with the 
standard mechanism for scarring is here discussed and analyzed. 
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1 Introduction 

The quantum calculation performed on the Bunimovitch 
stadium billiard [1] by McDonald and Kaufmann [2] showed, 
quite surprisingly, that some of the eigenfunctions present 
anomalously high density along certain (unstable) peri- 
odic orbits (PO) of the system. This example demon- 
strates in a very dramatic way the existence of a corre- 
spondence between quantum and classical mechanics [SJ. 

Five years later, Heller |4] established the basis for the 
corresponding theory [5] using semiclassical arguments, 
and coined the term "scarring" to name this phenomenon. 
In this seminal work, the physical mechanism for this en- 
hanced localization was explained, by considering the dy- 
namical evolution of a localized wave packet, |(/>(0)}, ini- 
tially launched in the vicinity of the scarring PO. The 
center of such packet follows for some time the classical 
path dictated by the PO, so that part of the initial quan- 
tum density distribution periodically returns next to the 
initial region. This fact originates recurrences in the self- 
correlation function C(t) = (0(O)|</>(f)), whose intensity 
decreases exponentially with time due to the competing 
mechanism induced by the unstable dynamics, character- 
ized by the Lyapunov exponent A, in the vicinity of the 
scarring PO. This discussion can be made more quantita- 
tive considering the relationship between the two associ- 
ated time scales, so that it can be stated that for a scar 
to be seen a necessary condition is given by the criterion 
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A/u << 1, where u> = 2tt/t is the frequency associated to 
the PO. 

Notice that the above arguments were formulated in 
the time domain, but they can be translated into the en- 
ergy domain by considering the spectrum associated to 
C{t) 

lEt ' n C(t). 



IiE) = h I 



dt i 



(1) 



This spectrum appears resolved in narrow peaks of width 
proportional to h\, occurring at intervals of hu>, as picto- 
rially depicted in Fig. 22 of Ref. [5]. Moreover, when this 
low resolution spectrum is projected in the system eigen- 
spectrum [6], the eigenstate intensities under these bands 
are enhanced by a factor of lu/ A with respect to a random 
distribution [7J. 

One important issue of this alternative view in the en- 
ergy domain is that it shows very clearly the relevance of 
the quantization of the action, S(E), along the scarring 
PO [3] in our problem. The energies, Ebs, for which this 
quantity fulfills the Bohr-Sommcrfcld (BS) quantization 
condition 

(2) 



S{E BS ) = 2Trh [n 



being n the corresponding quantum number and v the 
Maslov index associated to the orbit, mark the position of 
the centers of the low resolution bands appearing in {J). 
The reason for this is clear, since condition © guaran- 
tees the coherence of the probability recurrences along the 
PO, which is responsible for a constructive interference of 
the returning wave, causing the building up of probabil- 
ity density constituting the scar. This way of proceeding- 
is similar to the way in which the optical problems, such 
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as interference, are treated. There, the notion of optical 
path, analogous to our action, followed by the different 
waves all propagating at the same energy, is used. In this 
way, any reference to time can be avoided. 

Another important contribution to the theory of scars 
is due to Bogomolny [8], who derived an explicit expres- 
sion for the PO contributions to the smoothed quantum 
probability density over small ranges of space and energy 
(i.e. average over a large number of eigenfunctions). A 
corresponding theory for Wigner functions was developed 
by Berry [9]. Something interesting to remark about this 
work in connection with the present paper is that its au- 
thor noticed the existence of fringes as one moves away 
from the scarring PO, either in position or energy. Also, 
Prado and Keating focussed recently on the influence of 
bifurcations on scarring [10j . something of outmost rele- 
vance for (realistic) mixed systems. 

On the experimental side, many practical applications 
in areas of much interest have been described in the liter- 
ature concerning scar theory. Among them, nanodevices 
[TT] . optical microcavities [Hj, or optical fibers [113] de- 
serve a special mention. 

As discussed above, reference [4] mainly focussed on 
the probability density recurrently returning next to the 
initial point at time intervals corresponding to the scarring 
PO period, paying no attention to the fate of the remain- 
ing density leaving the PO dynamically unstable region. 
Obviously, this density can be thought as formed by differ- 
ent pieces, each of which would ramble through the avail- 
able phase space by a different pathway, each one of them 
with its own characteristic complexity. Actually, by sum- 
ming up all these contributions Heller and Tomsovic p3] 
were able to construct a semiclassical approximation to 
wave functions valid past the limit of the Ehrenfest time 
|15j , something that came as a surprise since at that point 
the underlying chaos has had time to develop structure at 
a scale much finer than h. 

The purpose of this paper is to study what happens 
with the portion of the quantum probability density first 
expelled from the vicinity of the scarring PO, but return- 
ing back to the initial region by the shortest circuit that 
is different from that defined by the scarring PO itself. In 
this sense, our goal is to present a complementary view to 
that provided by the standard scarring mechanism. The 
phenomenon that will be described constitutes a second 
order effect, since the amount of quantum density involved 
is orders of magnitude smaller than that transported along 
the scarring PO. The key point is to consider the prob- 
lem in the energy domain, so that the importance of ho- 
moclinic excursions related to the scarring PO (with the 
associated propagation of quantum waves) can be quan- 
titatively analyzed, and the relevance of the quantization 
of these phase space circuits gauged. 

The organization of this paper is as follows. In Sect. [5] 
we briefly described the model used in our calculations. In 
Sect. [3] we discuss some aspects of the classical dynamics 
of this system which are relevant to our work. We con- 
tinue with Sect. [4l in which we present and analyzed our 
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Fig. 1. Desymmetrized Bunimovitch stadium billiard and def- 
inition of the Birkhoff coordinates (q,p) used in our study. We 
focus our attention on the dynamics influenced by the horizon- 
tal periodic orbit drawn in dashed line. 

quantum results. We summarize our main conclusions in 
Sect. [5J 



2 Model 

The model used in our calculations consist of a particle of 
mass M — Ti 2 /2 confined in a desymmetrized Bunimovitch 
stadium billiard, defined by the radius of the circular part, 
r = 1, and the enclosed area, 1 + 7r/4, as shown in Fig. [1] 
Note that in this units system the energy of the particle is 
simply given by the square of the wave number, E = k 2 . 
This billiard constitutes a standard paradigm in the field 
of quantum chaos, and it has been shown to display all 
characteristics (ergodic, mixing, . . . ) of the chaotic hier- 
archy at classical level [IB] . 

To study the classical and quantum dynamics of this 
system, Birkhoff coordinates (on the boundary) [16] . (q,p), 
are customarily used. They can be seen superimposed in 
Fig. [T] Coordinate q, is the arc length coordinate, usu- 
ally measured starting from the upper left corner of the 
stadium, and p = p • t/|p| is the fraction of tangential 
momentum. 



3 Classical dynamical aspects 

In this paper, we will focus our attention on the dynamics 
influenced by the horizontal PO, which coincides in our 
case with the axis for the x coordinate. See dashed line in 
Fig.ffl 

A suitable phase space picture of the associated dy- 
namics can be obtained by defining a Poincare surface of 
section consisting on plotting the coordinates (q,p) when 
the trajectory hits the boundary (axis excluded to recon- 
cile with the dynamics in the full version of the stadium) . 
This is presented in Fig. [51 where some elements relevant 
to our work are shown. For example, the horizontal PO 
renders a fixed point located at (q,p) = (1 + n/2, 0), that 
appears labelled (P) in the plot. From this point two man- 
ifolds, one unstable and the other stable, emanate. They 
are plotted, with full and dashed line respectively, in the 
figure. As can be seen, these two manifolds first cross at 
the two primary homoclinic points, labelled hoi and I102 
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Fig. 2. (Left) Detail of the phase space portrait showing the 
relevant elements related with the dynamics of the horizontal 
periodic orbit shown in Fig. [T] 

Label (P) indicates the position of the corresponding fixed 
point. The associated unstable and stable manifolds are rep- 
resented in full and dashed line, respectively. They cross, with 
different topology, at the primary homoclinic points hoi and 
ho2, which map into hoi, ho", . . ., and ho 2 , ho 2 ', . . ., respec- 
tively. 

(Center and right) Primary homoclinic areas Shoi and Sho 2 
associated to the homoclinic circuit for the horizontal periodic 
orbit. 



respectively in the figure, and then at hoi,ho^,... and 
ho2, ho 2 , . . . (and their reflections with respect to the line 
p = 0). These crossings form lobes (the first of which 
is the only one shown in the figure) that by subsequent 
propagations with the Poincare map form the well known 
homoclinic tangle. This tangle organizes all the dynamical 
complexity of the system at the classical level, as shown in 
the early work of Poincare. These two infinite sequences 
of crossing points define the so-called primary homoclinic 
orbits, two in this case. As can be seen, they define in 
phase space two areas, 5h Q1 and >!5ho 2 , which are shown 
as shaded regions in the central and rightmost panels of 
Fig. [2l respectively. 

It is also interesting to examine these orbits in config- 
uration space. Since these orbits corresponds to infinitely 
long trajectories, accumulating at both ends in the hori- 
zontal PO, a complete calculation of such orbits it is not 
feasible. Accordingly, we will present here only a short- 
ened version of them that captures, however, the main 
dynamical part of it. These surrogates of the full orbits 
are computed in our case in the following way. We start 
by locating very accurately the position in phase space of 
the first two crossings of the unstable and stable mani- 
folds emanating from the fixed point (P); this giving the 
position of the two primary homoclinic points hoi and 
ho2- A piece of the required orbits is then calculated by 
propagating forward in time the corresponding trajecto- 
ries until they cross several more times (6 in our case) the 
SOS. The corresponding results are shown in Fig. [3] Obvi- 
ously, longer versions of the whole orbits can be produced 
by simply enlarging the above described procedure. This 



would involve taking the initial point closer to the fixed 
point and ending the trajectory at a later time, so that 
more points on the SOS at both ends of our previous ap- 
proximation to the full orbit are collected. However, this 
would only complicate visually our figures with more cum- 
bersome orbits, adding very little to the dynamics that we 
are now analyzing. 

4 Quantum mechanical aspects 

All quantum calculations performed in this work are car- 
ried out by imposing Dirichlet conditions on the stadium 
boundary and Neumman conditions on the (horizontal 
and vertical) axes. 

To perform our study we first need a method (the tool) 
to systematically construct a series of scar wave functions 
highly localized on a given PO. Let us remark in this re- 
spect that the appearance of scarred eigenfunctions be- 
comes less frequent as the energy increases, being van- 
ishingly small in the semiclassical limit. Accordingly, the 
scar wave functions that we need are, in general, non- 
stationary and can be obtained, for example, as linear 
combination of the eigenstates of the system, as described 
by Bogomolny [8]. 

Moreover, very sophisticated methods for performing 
this task has also been reported in the literature [6l [L7llT8j . 
In this work we use the method developed in Ref. |18j . 
which can be summarized as follows. We start from "tube 
wave functions", |</>t u be), covering the configuration region 
around the PO. They are obtained from the semiclassical 
theory of resonances developed by Vergini and Carlo (TTJ 
H"9] . and are calculated at the BS quantized energies [see 
eq. (|2])], which for our model can be expressed, in terms 
of the wave numbers, as 



with Lh = 4 (length of the horizontal PO) and v = 3. 
After that, we define a more elaborated scar function by 
incorporating the associated averaged [5] dynamical infor- 
mation up to a given time, T, and in the energy window 
Ebs ih/T. This is done by using the expression 

l^car) = dt COS g*(-^ BS —H)t/h (4) 

where the cosine filter is introduced to minimize the en- 
ergy dispersion of the resulting functions. 




Fig. 3. Primary homoclinic orbits passing respectively by hoi 
(a), and ho2 (b). The position in configuration space of these 
two phase space points has also been indicated in the plots. 
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Fig. 4. Scar function localized along the horizontal periodic or- 
bit corresponding to the quantum number uh = 54, in config- 
uration (a)-(b) and phase space (c)-(d). It has been computed 
with the aid of Eq. ([2]) using two different values of the time 
parameter T: (a)-(c) ts, and (c) 1.5tE, being tE the Ehrenfest 
time 1151. 



To monitor the different properties and characteristics 
of our scar functions, specially those aspects concerning 
the dynamical information carried by them, we resort to 
quantum surfaces of section (QSOS) obtained from suit- 
able (quasi)probability density distributions, in our case 
the Husimi function [20] . Phase space representations for 
the wave functions described in this section can be de- 
fined in a number of ways. In this paper we will follow the 
procedure described in Ref. [21], which relies on the use 
of normal derivatives of the wave functions evaluated at 
specific points on the boundary of the billiard, (q,p). The 
coherent states necessary to construct such representation 
(Husimi functions) are then defined as: 



1 



1/4 



exp 



1 

2^2 



(q - qf + ip(q - q) 



This expression corresponds to a boundary wave packet 
centered at the point (<?,p) on the surface of section, and 
represents a bounce off a given boundary point with a 
specific tangential momentum. Then, for a wave function 
with normal derivative on the boundary (j){q) (extended 
periodically to the real line), the corresponding Husimi 
function is given by 



H(q,p) = |(G 9 ,p|^ S car)r 



(6) 



As an example, we present in Fig. 01 both in config- 
uration (a)-(b) and phase space (c)-(d), the scar wave 
function associated to the wave number k = 86.01, value 
that been obtained using the BS quantization condition 
([3]) with nn — 54. Two values of the cut-off parameter, 
T, appearing in Eq. l[4"]). have been used, being those equal 
to T = tE for (a)-(c) and to 1.5i£ for (b)-(d), respec- 
tively. As can be seen in parts (a) and (b), the scar wave 
functions appear highly localized along the horizontal PO. 
Moreover, it should be remarked that the "quality" of 
these wave functions is very high since, for example, the 
focal point [5] existing on the horizontal axis at approxi- 
mately 3/4 of the PO total length is clearly visible in both 
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Fig. 5. (color online) Logarithm of the amplitude of the 
Husimi based quantum surface of section for the scar wave 
functions evaluated at points P (red full line), hoi (green 
dashed line) and ho2 (blue dashed-dot line), as a function of 
the corresponding wave number, k. 



plots. The reason for this has been thoroughly discussed 
in Ref. [T7], and comes from the fact that the dynamics 
in the vicinity of the scarring PO has been included in 
eq. ([!]). This is also illustrated in the results on the other 
two panels, (c)-(d), of Fig. 0] As can be seen there, the 
corresponding Husimi based QSOS is not only localized on 
the fixed point (P), corresponding to the horizontal PO, 
but also spreads in a substantial way along the paths of 
the associated manifolds. Moreover, this function spreads 
further along these manifolds as the time parameter T is 
increased. Actually, for the lowest value considered here, 
T = tE, the probability distribution only covers the linear 
part of the manifolds, while for the larger value T = \Me 
it extends much further, getting even to the region where 
the primary homoclinic points, hoi and ho2, are located. 

Let us now return to the main point addressed in this 
paper, namely, what it is the fate of the probability density 
which after being launched from a point in the vicinity of 
the horizontal PO, instead of returning next to it in the 
first recurrences, it is thrown out of this region. As will 
be shown, this can be unveiled by analyzing the behav- 
ior of the probability density for different values of the 
quantum number nn, or equivalently of the correspond- 
ing BS quantized wave number, k. Instead of considering 
the whole distribution, we will simplify our calculation by 
focussing the attention into a single relevant point. Fixed 
point (P) seems to be the most representative one, since 
it is in its vicinity that the density returning by the re- 
current mechanism discussed by Heller in his 1984 paper 
4J, is coherently collected. Let us remark once again that 
this is so because in the way we perform our calculations, 
the action along the PO circuit is quantized. 

The results for 91 scar wave functions corresponding 
to n = 10 - 100 (k = 16.88 - 158.25) are shown in red full 
line in Fig. [5j In it we see that, superimposed to a slightly 
increasing overall tendency, there is a small oscillatory be- 
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Fig. 6. (color online) Fourier transform of the amplitudes in 
phase space presented in Fig.[S]after having removed the slowly 
varying contribution as described in the text. 

havior. Although this is not well visible in the scale of the 
figure, it is obvious when one considers its Fourier trans- 
form, which is shown in Fig. [6] At this respect, notice 
that the conjugate variable used, S, has units of length, 
and plays the role of an action in our calculations due to 
the unit system that we have chosen. 

As can be seen, several peaks contribute to the oscil- 
latory behavior of logH computed at phase space point 
(P), being those at S = 0.633 and 1.007 the most impor- 
tant ones. The origin of this regularity can be now quite 
clearly understood. In the first place, the overall value in 
the QSOS function, H, comes from the scarring mech- 
anism by recurrences along the (horizontal in this case) 
PO, as described by Heller in Ref. [4]. The actual value 
may change slightly with n#, but this variation is small, 
since the essential constructive mechanism induced by the 
PO recurrences is the same in all cases. However, and as 
proposed by us, we have a second, superimposed mecha- 
nism. This consists in recurrences due to longer, and in- 
creasingly more complicated paths in phase space. From 
all of them, the shortest possible are those taking place 
along the homoclinic circuits defined by the two primary 
homoclinic orbits described in Sect. [3] This second inter- 
ference mechanism is controlled by the value of a different 
action, and then the corresponding quantization condition 
does not coincide, in general, with that of the scarring PO. 
However, every so often the two quantization conditions 
will take place at very close energies; this will correspond 
to a maximum in the value of the Husimi function, H , 
computed at the phase point (P). Obviously, at the mid- 
point between two consecutive quantized energies of these, 
we will be in an "antiquantization" point, where the inter- 
ference of the quantum density when arriving at (P), after 
a recurrence along the homoclinic circuit, will be destruc- 
tive, and then we will find a minimum in H . Intermediate 
situation will then be found at the energies in the range 



defined by the two quantized energies, and this is the ori- 
gin of the oscillatory behavior we are observing in the red 
full line of Fig. [5] 

To confirm numerically this interpretation, and explain 
why we are observing two principal frequencies instead of 
one, let us consider now the quantization condition(s) of 
the homoclinic circuit denned by the manifolds emanating 
from (P) (see two rightmost panels in Fig. [2]). As men- 
tioned before the homoclinic orbits denning these circuits 
are infinite in length, something which makes difficult to 
derive the quantization condition we are looking for. The 
way to deal with this problem was discussed in detail in 
Rcf . 22] , and it is based on an idea of Ozorio de Almeida 
[23] , when he tried to quantized the so-called "homoclinic 
torus" . Basically, it consists of taking advantage of the 
fact that the homoclinic orbits both leave and goes back 
asymptotically to the horizontal orbit (P). In this way, 
instead of calculating the action along the homoclinic or- 
bits, one can just compute the difference with the hori- 
zontal one, which is already quantized. Accordingly (see 
full details in [22]) 

7T 

kS h 0i - ^ ^ho, = 27T71, 1 = 1,2, (7) 

which is valid for both primary homoclinic orbits/circuits, 
hoi and I102. In our case it is not difficult to show that 
i\ 0l = — 1 and fj 102 = u 2J], and the actual values of 5h 0l 
and Sho 2 can be numerically calculated with the aid of 
Fig.H resulting equal to -3.368 390 45 and -2.991 142 22, 
respectively. These two quantities are in excellent agree- 
ment with the two values, 0.633 and 1.007, found in the 
Fourier transform of H evaluated at (P) (see Fig. [6]), after 
subtracting the length of the horizontal orbit, Lh = 4. 

To conclude this section it should be pointed out that 
similar results are obtained if apply the same type of anal- 
ysis presented above to the fluctuations of the Husimi 
function evaluated at the two primary homoclinic points, 
hoi and I102 . The results are shown in Fig. [5] with green 
dashed and blue dashed-dot lines, respectively. In these 
two cases, the oscillations are much more pronounced, 
but the Fourier transform (Fig. [5]) shows identical results. 
Namely, for the first curve a single big peak, correspond- 
ing to S = 0.633 exists, while the oscillations for the sec- 
ond one are fully accounted for with a single frequency 
at 1.007. Both results are in perfect agreement with the 
results and conclusions discussed before. 



5 Conclusions 

In this paper, we have revisited the scarring mechanism 
proposed by Heller in his 1984 paper [4], in which he ana- 
lyzed the role played by recurrences along classical unsta- 
ble POs. We perform our study in the energy rather than 
in the time domain, as it is customary when arguments 
based on the self-correlation function, C(t), are used. For 
this purpose, sophisticated scarred wave functions, highly 
localized along POs, have been used as a tool for our anal- 
ysis. Special attention is paid to the fate of the quantum 
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probability density first expelled from the unstable vicin- 
ity of the scarring PO, showing the relevance of the as- 
sociated homoclinic circuits for the coherent, i.e. causing 
constructive interference, transport of probability density 
in classically chaotic systems. 
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